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Abstract 

We consider the hypothesis that neutrino oscillation data can be 
explained if the gravitational couplings of (massless or degenerate 
mass) neutrinos are flavour non-diagonal, in violation of the equiv- 
alence principle. We analyze the various neutrino oscillation labo- 
ratory experimental data including the recent LSND observations to 
constrain the relevant parameter space. We find that there is no al- 
lowed region of parameter space which can explain the existing data, 
implying that the LSND result cannot be explained by oscillations of 
degenerate-mass neutrinos due to equivalence principle violations. 
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Empirical evidence supporting neutrino-flavour oscillations continues to 
mount 0. At present there are four different solar neutrino experiments ||, 
each using distinct detection techniques, that consistently find a discrepancy 
between the measured solar v e flux and that predicted by solar models [||]. 
There are also a number of experiments on atmospheric neutrinos |3J which 
find that the ratio of the flux of to v e is significantly smaller than one 
would expect from standard particle physics models [f|. Most recently, the 
Liquid Scintillator Neutrino Detector (LSND) group has recently announced 
an excess VI events between the energy range 36 and 60 MeV || . If this excess 
is due to the z^J — >■ Z^T oscillation, then it implies an oscillation probability of 
(0.34l ;i 8 ± 0.07)%. The distance traversed by the V^ before being detected 
as a Pe is about 30 metres. 

Mechanisms underlying neutrino oscillation typically assume that neu- 
trinos have nondegenerate masses, following the original suggestion by Pon- 
tecorvo 0. In this scenario the weak interaction eigenstates of neutrinos are 
distinct from their mass eigenstates, thereby permitting oscillations between 
the various flavours. 

An alternative neutrino oscillation mechanism was proposed more re- 
cently by Gasperini || (and independently by Halprin and Leung ||), in 
which neutrino weak interaction eigenstates are distinct from their gravita- 
tional eigenstates. This mechanism (later referred to as the VEP mechanism 
[ PHI ) does not require neutrinos to have nonzero masses; instead neutrino 
oscillations occur in this mechanism due to an assumed flavour non-diagonal 
coupling of neutrinos to gravity, in violation of the equivalence principle. 

From this viewpoint, neutrino oscillation experiments furnish us with a 
test of the equivalence principle. The VEP mechanism has been explored in 
a number of papers [|ll|] as a possible explanation of solar neutrino data. A 
recent analysis jTIJ has shown that, in the context of a two-flavour model, 
there are small allowed regions of parameter space at both small and large 
VEP mixing angles which are compatible with present day solar observations. 
Extension to a full three-flavour model indicates that the allowed regions of 
parameter space can widen due to mixing with a third flavour [|T2|]. 

In the present paper we consider the VEP mechanism in the context of 
laboratory searches for neutrino oscillations. We first show that LSND data 
itself can be explained by neutrinos of degenerate or zero mass with flavour 
non-diagonal gravitational couplings. We then carry our analysis further 
to include other laboratory experiments [13, |T^, I!J, which also constrain 



1 



the allowed VEP parameter space. We find that the combination of these 
constraints rules out any violation of the equivalence principle in the e — /i 
sector, implying that gravity couples to u e the same way as v^. 

As a consequence, in the absence of other physical mechanisms (such as 
lepton number violation), previous accelerator data in conjunction with the 
LSND results can only be explained by assuming neutrinos have differing 
masses. A repetition of the above analysis in this case would then imply a 
new bound on the parameter space of the VEP mechanism ||. 

In this article we shall not argue in what circumstances the equivalence 
principle might be violated. Rather, we take a phenomenological approach to 
this problem and try to constrain the parameter space only from an analysis 
of the existing data on neutrino oscillations. For the sake of simplicity we 
work in the two generation scenario. At the end we shall comment on the 
role other oscillation data plays in constraining VEP. 

We turn now to the question of whether or not the LSND results can 
be understood solely in the context of the VEP mechanism. In this mecha- 
nism, the gravitational eigenstates \vq) = (vig, v 2g) are related to the weak 
eigenstates \vw) = (i/ e , i^) by an SO (2) rotation R(9a) 

Ww) = R(0g)Wg) (1) 

where Qq is the mixing angle. The gravitational eigenstates are solutions to 
the Dirac equation in a Schwarzchild background. For a spherically symmet- 
ric metric, choosing the trajectory of the neutrino in the radial direction, we 
can write down the diagonal components of the Hamiltonian, 

fl? = -2|#r)|3i(l + / 4 ) (2) 

which governs the evolution of the neutrinos. Here the fi are the flavour 
dependent gravitational parameters, which determine the magnitude of the 
violation of the Weak Equivalence Principle. The evolution of the weak 
eigenstates will be governed by the equation 

dr \p^J lY( n V§sm20 G cos2# G J \u^J v ; 

where \4>{r)\ is the Newtonian gravitational potential and Af = f 2 — f\. If 
the equivalence principle is not violated then fx — / 2 . 
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In this paper we shall be discussing small scale terrestrial laboratory ex- 
periments, for which <p(r) may taken to be constant. Although a natural 
choice for <f> would be the earth's gravitational potential (~ O(10 -9 )), an- 
other choice is to consider the potential due to all forms of distant matter. 
The dominant contribution is from the local supercluster which has been es- 



timated to be 3 x 10" 5 17 . For our purposes the choice of <fi is irrelevant, 
since to find the allowed parameter space we shall consider \(j)\Af as the rel- 
evant parameter. Particular limits on Af that arise from given experiments 
may be found by substituting the above values for (/). 

Consider a beam of muon neutrinos that traverses a distance of L me- 
ters. The probability that a muon neutrino will get converted to an electron 
neutrino is given by 

P(v„ -> v e ) = sin 2 29 G sin 2 — , (4) 

where X G = E\</>(r)\Af * s ^ ne oscillation length. Although this oscillation mech- 
anism has a number of similarities with that for neutrinos of non-degenerate 
mass, the oscillation length has a markedly different energy dependence, 
varying inversely with energy in (^), whereas it is proportional to energy in 
the massive case [0. As we shall see, although there exists a fairly wide pa- 
rameter space for neutrinos of non-degenerate mass which can explain all the 
laboratory experiments, the VEP mechanism for degenerate mass neutrinos 
cannot explain all these experiments. 

To analyse any neutrino oscillation data, it is useful to divide the param- 
eter space into three regions: L << X, L » X and L ~ A. where A is the 
oscillation length. In the large A regime, sin 2 29 can be as large as unity, and 
the data constrain A as a function of the mixing angle 9. As 9 decreases, 
the maximum excursion for this part of the curve occurs at approximately 
sin 2 29 = (P), for which A = 2L. In the case of the LSND experiment, if 
one interprets the excess TT e events as neutrino oscillations, then the above 
oscillation probability would mean both lower and upper bounds on X G , and 
hence on \(p\Af. In the small A region large numbers of oscillations take 
place before the beam reaches the detector, and so the value of sin 2 ^ may 
be assumed to be its average value |. In this regime, observation of a small 
(or null) neutrino oscillation probability puts both lower and upper (or just 
upper) bounds on sin 2 29 for large \<p\Af. In the intermediate region the 
oscillation length is comparable to the travel length L and so sin 2 ^ varies 
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slowly between and 1, putting limits on sin 2 26 that are quite sensitive to 
the actual value of \(fi\Af. 

For the LSND experiment, the oscillation probability (^) in the VEP 
mechanism can be simplified to 

P = sin 2 2# G sin 2 (7.62 x lO 15 |0(r)|A/) (5) 

where we have taken (L) = 30 meters and the average neutrino energy is 
taken to be 50 MeV. The change on the bounds on \<j)(r)\Af for large sin 2 26 q 
due to the uncertainty in the distance or neutrino energy are negligible. For 
sin 2 2#G = 1, the allowed region (95% C.L. from LSND |J) for the violation 
of the equivalence principle is 

9.76 x 10" 18 > \<f>(r)\Af > 5.03 x 10~ 18 . (6) 

The small A region occurs when 

\(j)\Af > 1.35 x 10~ 16 (7) 

for which the LSND data yields the bound 

.0029 < sin 2 26 G < .011. (8) 

Thus with the above bounds on the VEP parameter \(p\Af, one can explain 
the LSND result. This would apparently mean that LSND result does not 
imply a non-trivial neutrino mass matrix. We shall now demonstrate that 
these allowed regions have already been ruled out by other accelerator ex- 
periments. 

The E776 experiment at the Brookhaven National Laboratory did not 
see any statistically significant excess number of v e (pg) events over the back- 
ground at a distance of 1 km from the source in a wide-band (z^j) beam. 
Most of the events are above 1 GeV and peaked around 1.4 GeV. Only 19 
events with an expected background of 25 ± 5 ± 3 ± 3 were observed and from 
this an upper limit on the probability of neutrino oscillations was determined 



at 90% C.L. (TJ 



For large sin 2 6q this gives a lower limit on the mass-squared difference 
consistent with the LSND result for a neutrino oscillation scenario of neutri- 
nos of non-degenerate mass [0] . However, because of the difference in energy 
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between the LSND and the E776 experiments, the 1 GeV null result of the 
E776 gets translated to 

\<f)\Af < 3.0 x 1(T 21 (9) 

in the VEP mechanism. This bound is not consistent with the LSND result 
(ID). In table 1, we present the upper and lower bounds on \<p\Af that are 
permitted within the limits of error from the E776 and LSND experiments 
respectively for several values of sin 2 26 g- Over the entire range we find that 
the upper bound on the \4>\Af from the E776 experiment is much less than 
the lower bound of |0|A/ as allowed by the LSND result. As a result there 
does not exist any region of the parameter space of the VEP mechanism 
which can explain both the LSND result as well as the E776 experiment for 
sin 2 26 G > .003. 

We now concentrate on the other regions of the parameter space. The 
small \q region (for which \<p\Af is so large that many oscillations occur 
within the experimental beam length) sets in for the E776 experiment at 

|0|A/ > 8.03 x 10" 20 (10) 

which is consistent with (|7|), and implies a bound of 

sin 2 2# G <.003 (11) 

on the mixing angle. Thus combining the allowed regions of the E776 exper- 
iment with the LSND result implies that both these experiments can only 
be explained by degenerate mass neutrinos through the VEP mechanism for 
very large \ fi\Af > 1.35 x 10" 16 and for 0.0029 < sin 2 29 G < 0.003. 

Even this marginally consistent result may be ruled out as follows. If we 
further argue that gravity couples with matter and antimatter in the same 
way (so that gravitational interactions conserve total lepton number), then 
we can compare the allowed parameter space of the VEP mechanism with the 

- * v e oscillation limit as obtained by the SKAT experiment at Serpukhov 



14], which provides the most stringent upper limit on the mixing angle in 



the small A region at 90% C.L. | 15| . SKAT measures the ratio of v e to 
induced charged current reactions as observed in a bubble chamber exposed 
to a wide band neutrino beam with energies between 3 and 30 GeV and a 
neutrino beam length of 270 meters. The small \q region sets in for the 
SKAT average beam energy at |0|A/ > 2.3 x 10~ 19 and the bound on the 
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Tabic 1: Bounds on |0(r)|A/ from E776 (upper) and LSND (lower) 





|0(r)|A/ 


sin 2 26 G 








upper bound from 


lower bound from 




E776 experiment 


LSND experiment 


0.002 


8.03E-20 


1.35E-16 


0.003 


6.06E-20 


1.02E-16 


0.004 


5.09E-20 


8.55E-17 


0.005 


4.48E-20 


7.52E-17 


0.006 


4.05E-20 


6.8E-17 


0.007 


3.72E-20 


6.25E-17 


0.008 


3.46E-20 


5.81E-17 


0.009 


3.25E-20 


5.46E-17 


0.01 


3.08E-20 


5.16E-17 


0.02 


2.15E-20 


3.6E-17 


0.03 


1.75E-20 


2.93E-17 


0.04 


1.51E-20 


2.53E-17 


0.05 


1.35E-20 


2.26E-17 


0.06 


1.23E-20 


2.06E-17 


0.07 


1.14E-20 


1.91E-17 


0.08 


1.06E-20 


1.78E-17 


0.09 


1.00E-20 


1.68E-17 


0.1 


9.51E-21 


1.6E-17 


0.2 


6.71E-21 


1.13E-17 


0.3 


5.48E-21 


9.19E-18 


0.4 


4.75E-21 


7.96E-18 


0.5 


4.24E-21 


7.12E-18 


0.6 


3.87E-21 


6.5E-18 


0.7 


3.59E-21 


6.02E-18 


0.8 


3.35E-21 


5.63E-18 


0.9 


3.16E-21 


5.31E-18 


1 


3.00E-21 


5.03E-18 
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mixing angle is sin 2 26g < .0025. This fully rules out the the VEP mechanism 
for degenerate mass neutrinos. 

We can thus conclude that the LSND result cannot be explained by neu- 
trinos of degenerate mass if other laboratory bounds on neutrino oscillations 
are taken into consideration. This result holds regardless of the gravitational 
potential at the earth's surface. In particular, for massless neutrinos there is 
no allowed region of the parameter space of the VEP mechanism which can 
explain all the experiments. 

If we admit the possibility of a non-trivial neutrino mass matrix, then 
all of the above experiments only put bounds on the VEP parameters. As 
discussed in ref . || , in the expression for the neutrino oscillation probability 
in a two-flavour scenario the oscillation length A is now a function of \4>\Af, 
Am 2 and two mixing angles, and the LSND experiment will limit a combi- 
nation of these parameters in the large oscillation region via an appropriate 
generalization of eqs. (|^) and (||). The other experiments will provide fur- 
ther constraints on the parameter space, and there will be some minimal 
Am 2 which is consistent with all empirical constraints. We intend to relate 
details of this analysis in a forthcoming paper. 

In summary, we have shown that the LSND result in conjunction with 
other laboratory experiments rules out the possibility of the VEP mechanism 
for neutrinos of degenerate mass (this includes massless neutrinos as well). 
This situation arises because of the particular energy dependence of neutrino 
oscillations in the VEP mechanism in combination with the differing neu- 
trino energies employed in each experiment. In the absence of other physical 
mechanisms for introconversion of neutrino species, these results imply that 
neutrinos must have different nonvanishing masses. A naturalness argument 
would then imply that if the gravitational couplings of v e and are the 
same, then the gravitational coupling of v T should be equal to these, making 
the VEP mechanism an unlikely candidate for a neutrino oscillation mech- 
anism. Once we admit the possibility that neutrinos are massive, there is 
little motivation to consider the VEP mechanism as the mechanism chiefly 
responsible for neutrino oscillations. Present experiments can at best put 
bounds on the relevant parameter space. 
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